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Abstract 

In 1676 Olaus R0mer presented the first observational evidence for a finite liglit velocity Cem- 
He formed his estimate by attributing the periodicahy varying discrepancy between the observed 
and expected occultation times of the Galilean satellite lo by its planetary host Jupiter to the 
time it takes light to cross Earth's orbital diameter. Given a stable celestial clock that can be 
observed in gravitational waves the same principle can be used to measure the propagation speed 
Cgw of gravitational radiation. Space-based "LISA" -like detectors will, and terrestrial LIGO-like 
detectors may, observe such clocks and thus be capable of directly measuring the propagation 
velocity of gravitational waves. In the case of space-based detectors the clocks will be galactic 
close white dwarf binary systems; in the case of terrestrial detectors, the most likely candidate 
clock is the periodic gravitational radiation from a rapidly rotating non-axisymmetric neutron 
star. Here we evaluate the accuracy that may be expected of such a R0mer-type measurement 
of Cgw by foreseeable future space-based and terrestrial detectors. For space-based, LISA-like 
detectors, periodic sources are plentiful: by the end of the first year of scientific operations a LISA- 
like detector will have measured Cgw to better than a part in a thousand. Periodic sources may 
not be accessible in terrestrial detectors available to us in the foreseeable future; however, if such 
a source is detected then with a year of observations we could measure Cgw to better than a part 
in a million. 
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I. INTRODUCTION 

Over the course of a Jovian synodic year the distance hght must transit between Earth 
and Jupiter varies by two astronomical units: approximately 3 x 10^ km. If we neglect the 
time required for light to cross Earth's orbit, the interval between events that are periodic 
at Jupiter will at Earth be observed at times that may vary from the expected by as much 
as 10^ s. It was this observed variation between the observed and expected occultations of 
the Galilean satellite lo that led Olaus R0mer to conclude that light has a finite propagation 
speed and to the first real measurement of light's velocity Cem [III2]- Galactic compact white 
dwarf binary systems, or rapidly rotating non-axisymmetric neutron stars, are similarly 
regular clocks whose periodic signal is transmitted to Earth via gravitational radiation. In 
the same way that R0mer was able to use observations of the discrepancies in the optically 
measured times of lo's occultations by Jupiter to measure the speed of light, so we may use 
gravitational-wave observations by space-based LISA-like detectors [3H7] of compact white 
dwarf binary systems, or by terrestrial LIGO-like detectors |5HTT] of rapidly rotating neutron 
stars, to measure the propagation speed Cg„ of gravitational waves. 

In general relativity theory gravitational waves and light wave both propagate on space- 
time null geodesies; correspondingly, there is no difference in their respective (vacuum) 
propagation speeds. A direct measurement of the gravitational-wave propagation speed is, 
thus, a "go/no-go" test of the theory. 

Measurement of the R0mer delay directly and unambiguously access the wave propaga- 
tion speed across Earth's orbital baseline. This stands in contrast to other proposed tests of 
general relativity whose results are sometimes discussed in terms of the gravitational-wave 
propagation speed, but whose interpretation in this way requires a theoretical model or phe- 
nomenological framework to relate the observation to the propagation speed. For example. 
Will [12j suggested searching for an anomalous (compared to general relativity's prediction) 
compression of the gravitational- wave signal from an inspiraling binary system. Such a 
compression could be interpreted as a frequency-dependent gravitational- wave propagation 
speed resulting from a non-zero "graviton-mass" . In a similar vein, Larson and Hiscock 
[13] . Cutler et al. [H] proposed measuring the phase difference between the binary's orbital 
phase at some fiducial time as determined optically and by gravitational-wave observations. 
The phase difference, relative to that predicted by general relativity, could then be inter- 



preted as differences in the gravitational-wave propagation speed at a frequency twice the 
binary orbital frequency. What is relevant is that, unlike the measurement described in this 
work, none of these other measurements directly accesses the gravitational-wave propagation 
speed: i.e., their interpretation in terms of the wave propagation speed requires a theory or 
phenomenological framework that relates the observed phenomena to the wave propagation 
speed. As shown by Carlip [15J in the context of a recent claim to have measured the prop- 
agation speed of the gravitational force [16] through its effect on light travel, change the 
theory or framework and the interpretation changes. 

In Section [IT] we estimate the precision to which we can measure the gravitational-wave 
propagation speed from multi-year observations of periodic gravitational waves. We assume 
here that the gravitational- wave frequency and sky location of the source are known a priori. 
For such cases, we show that the Fisher matrix estimate of the uncertainty can be expressed 
very simply in terms of the source parameters and orbital radius of the Earth's motion 
around the Sun, valid for all ground-based and proposed space-based detectors. Details 
specific to the detectors, such as antenna pattern functions, cancel out when the uncertainty 
is expressed in terms of the signal-to- noise ratio of the measurement. We also discuss the 
complications introduced if we relax the assumption of a priori knowledge of the source 
frequency and sky location of the gravitational-wave source. Although the calculation is 
more complicated for this case, the Fisher matrix formalism can still be used to estimate 
the uncertainty in the gravitational-wave propagation speed as a function of the source 
parameters and detector geometry. In Section III we use the general result of Section [IT] to 



obtain numerical values for "3cr" fractional uncertainties in Cg„ for observations in terrestrial 
and space-based detectors. We also summarize our conclusions. 



II. ESTIMATING THE MEASUREMENT PRECISION OF THE 
GRAVITATIONAL- WAVE PROPAGATION SPEED 

A. Introduction 

We use the Fisher Information Matrix formalism [T71 [TS] to estimate the precision with 
which we can estimate the gravitational-wave propagation speed from multi-year observa- 
tions of periodic gravitational waves whose frequency and propagation direction are known 



a priori. For example, in the case of observations in a ground-based detector or detector 
network, the source may be a rapidly rotating neutron star whose rotational frequency and 
sky location are known from observation of its radio pulses; or, in the case of observations 
made with a space-based detector, the source may be a close white dwarf binary system 
that has also been observed optically. 

For any monocrhomatic source we may write the TT gauge gravitational wave strain at 
time t and location x as 

h{t,x) = h^{t,x)e^ + hy^{t,x)ey^ , (la) 

where e+ and e^ are orthogonal polarization tensors, fixed in inertial space, and 

h+{t, x) = H+ cos [27r/g„M + $+] , (lb) 

/ix (t, x) = Hy, cos [27r/g„M + $x] , (Ic) 

u = t — ^ . (Id) 

Cem(l + e) 

Here e is the fractional difference between light speed (cem) and the gravitational-wave 
propagation speed (cgw = Cem(l + e)), k is the wave-propagation direction and iJ+, H^, $+ 
and $x are constants determined by the source orientation with respect to e^ and Cx- 

For detectors that are small compared to the observed radiation wavelength [23j we may 
write the detector response to the incident (/i+,/ix) as 

r{t) = F+{t)h+{t, x{t)) + Fx (t)/ix (t, x{t)) , (2) 

where x{t) is the detector location. The functions -F+ and Fx are determined by the projec- 
tion of the detector's antenna pattern on the wave polarization tensors. Both terrestrial and 
space-based LISA-like detectors are constantly changing their orientation with respect to e_|_ 
and Cx (in the case of terrestrial detectors owing to Earth's diurnal motion about its rota- 
tion axis, and in the case of space-based detectors owing to the science-craft constellation's 
annual motion about about Sol); correspondingly, F+ and Fx are time dependent. 

B. Known source location and frequency 

The Fisher Information matrix X has elements 



where T is the observation duration, o"„ is the detector noise power spectral density at the 
gravitational wave frequency /g„, and 9 denotes the parameter vector (e, if+, ifx, *&+) "^x)- 
The partial derivatives of r with respect our parameterization are 

— = -27r/g^ {F+H+ sin [27r/g^M + $+] + F^H^ sin [27r/g„M + $x]} , (4a) 

dr 

: F+ COS (27r/g„M + $+) , (4b) 



dH+ 
dr 

dr 

dr 



FxCos(27r/g„M + <l'x) , (4c) 

-F+H+ sin (27r/gwM + $+) , (4d) 

-Fx i/x sin (27r/g„M + $x ) • (4e) 

For all cases of interest the gravitational-wave detectors follow Earth's orbit about Sol; 
corresp ondingly, 

k ■ X = (-Rau cos 9) cos{ujQt — (f)) , (5) 

where i?au is Earth's orbital radius (1 au), Uq is the detector angular velocity in its orbital 
motion about Sol (27r/yr), 9 is the ecliptic latitude, and is the azimuthal angle of the source 
with respect to Earth's orbital position at t = 0. (See Figure [I]) (The small displacement 
d oi a terrestrial detector away from the Earth's orbital path about Sol introduces an order 
d/Ra.u ~ 0.25% correction, which we ignore.) 

To evaluate the components of the Fisher matrix we take advantage of the sinusoidal 
periodicity of F+, F^, /i+ and h^ and focus on observations that are integer multiples of a 
year duration. Noting that 

^0 < 27r/gw < Ccm/d , (6) 

W© < Cem/^au , (7) 

the integrals for the Fisher matrix elements X^j for T > 1 yr quickly simplify to 

_ (27r/g„i?auPcos6')2 

1,^ = forje{ff+,//>,,«+,*>,}. (9) 

where p^ is the (power) signal-to-noise ratio 

p' = ^,fr'{t)dt. (10) 

•^"n Jo 
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FIG. 1: The relevant geometric quantities used in the calculation: k is the unit vector pointing 
in the direction of wave propagation; is the ecliptic latitude (i.e., the angle that —k makes with 
the plane of the ecliptic); is the azimuthal angle of the source with respect to the Earth's orbital 
position at t = 0. The detector antenna pattern functions -F+ and -Fx from Eq. pi) are defined 
with respective to the polarization tensors constructed from u and v, which are proportional to the 
unit vectors (p and 9, respectively. 

Correspondingly, at the level of the Cramer-Rao bound there is no co-variance between the 
uncertainty in our estimate of e and any of the other problem parameters. The expected 
variance of the estimate for e is thus 
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(11) 
(12) 



This result is valid for observations in all ground-based detectors and all proposed space- 
based detectors, whether Earth- or solar-orbiting. It is also valid for detector arrays where 
the data are combined coherently as described in, e.g., [191 ED]- (In the case of detector 
arrays p^ is the array power signal-to-noise ratio.) Details specific to the detectors, such as 
the antenna pattern functions F+ and F-^ , cancel out when the uncertainty is expressed in 
terms of the signal-to-noise p. 



C. Unknown source location and/or frequency 

If the source frequency, sky location, or both are not known a priori one needs to enlarge 
the parameter vector 9 to include the additional unknowns: e.g., the frequency /g„ and/or 
the source location on the sky {6, (p). The Fisher matrix dimensionality thus expands to 
include terms involving partial derivatives 

flT' I n ' IT \ 

-27r i TTT^ {F+H+ sin [2'Kf^^u + $+] + F^H^ sin [27r/g„u + $x]} , 



<9/g„ V Cem(l + e 

dr '^nfg^ d{k ■ x 



(13) 



de Cemfl + e) 86 



{F+H+ sin [27r/g^M + $+] + F^H^ sin [27r/g„M + <l>x]} 



+ l^H+ cos [27r/g^M + $_^] + -^ Hy, cos [27r/g„M + $x] | , (14) 

{F+H+ sin [27r/g^n + $+] + F,,H,, sin [27r/g„M + $x]} 



dr 27r/gw d{k ■ x 



Cem(l + e 



where 



+ I^H+ cos [27r/g„M + $+] + ^ ff, cos [27r/g„M + $x] | , (15) 



= — -Rau sin 9 cos(u;0t — 0) , (16) 



9(A; ■ X 



-i?au COS ^ sin(ci;0t — 0) . (17) 



one can 



Comparing these expressions with the partial derivative dr/de from Equation 4a, 
see that the off-diagonal Fisher matrix elements I^j for j G {/gw, 9, 0} are non-zero. Cor- 
respondingly, the elements of the covariance matrix (X^^)^^ are no longer trivial and z/^ no 
longer simply expressed. How well we can estimate the propagation speed of gravitational 
waves using observations of periodic sources whose frequency or sky location are not known 
a priori is the subject of work in-progress. 

III. DISCUSSION 

As described here, to measure the gravitational-wave propagation speed from the R0mer 
delay it is necessary to monitor a periodic source of gravitational waves, whose position 
on the sky and radiation frequency is known, for a year or more. For terrestrial detectors 
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such a source might be a radio pulsar that also radiates gravitationally. For such sources, 



equation 11 may be written as 

„,.(3.2xlO-')^(Hli!)'(Hy(l + .,Sec^9; (18) 

i.e., the "3cr" fractional uncertainty in the measurement of the gravitational- wave propaga- 
tion speed arising from a one or more year observation of a 100 Hz gravitational- wave source 
situated on the ecliptic plane is 10~^(10/p). Since a signal-to-noise p ~ 10 is typically taken 
as the threshold for source detection in a ground-based detector or detector network, if a pe- 
riodic source is observed a measurement of Cgw to 3a precision, 300(10/p) m s~^, will follow 
shortly. 

There are no reliable predictions for the gravitational-wave amplitude associated with 
rapidly rotating neutron stars. Mass asymmetries — "mountains" — are limited in size 
by the tensile strength of the neutron star crust, while the potential for fluid circulation 
instabilities (r-modes) to lead to gravitational- wave emission depends on the (temperature 
dependent) neutron star surface "ocean" shear and bulk viscosities [2TI §7.3]. It may well 
be the case that neutron star crusts cannot support a sufficiently large asymmetry to be 
observable gravitational-wave sources, or that the neutron star fluid viscosity is always so 
great as to stabilize neutron star fluid r-modes. Likewise, it may be that circumstances 
can be contrived that lead neutron stars to be strong radiation sources for ground-based 
detectors, but that there is no natural mechanism for creating or placing the neutron star 
into such states. Thus, while a sensitive measurement of the gravitational-wave speed is 
possible with ground-based detectors, carrying it out depends on the observation of a type 
of source that may not be available to us. 

Strong sources of periodic gravitational waves, in the form of galactic white dwarf binary 
systems, are both certain and plentiful for any of the proposed "LISA" -like space-based 



gravitational wave detectors [3]. For such sources. Equation 11 is conveniently written as 

, 4,2 /10mHz\ VlOO\\ x4 9. 
u, = {3.2 X IQ-y i— j ( j {l + e)\ec^e. (19) 

An amplitude signal-to-noise of 100 in a one-year observation is the minimum expected for a 
typical "verification binary" in a space-based detector; correspondingly, the "3o"" fractional 
uncertainty in the measurement of the gravitational-wave propagation speed arising from a 
one or more year observation of a 10 mHz gravitational-wave source situated on the ecliptic 
plane is a quite respectable 10~^(100/p). 



At present, then, we find ourselves in an odd position. With the observation of a periodic 
gravitational- wave source, we know how to make a direct, accurate and unambiguous mea- 
sure of the gravitational- wave propagation speed and, from it, a "go/no-go" test of general 
relativity theory. On the one hand, for existing or foreseeable future ground-based detec- 
tors there is no guarantee that there will, or — indeed — can, exist any source that will 
enable the measurement. On the other hand, there are scores sources, already identified, 
that are accessible to a space-based detector that would enable such a measurement but, 
despite the strong recommendation of the United States National Research Council |22], 
NASA abandoned its committment to the decade-long ESA/NASA partnership that would 
have led to the construction of such an observatory and no such project is currently planned 
by either agency. We can only hope that the feasibility of accurately and unambiguously 
testing general relativity — by means such as described here — will strengthen the case for 
reviving a LISA-like gravitational wave observatory in the near future. 
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